A PROOF OF THE MAIN THEOREM ON BEZOUTIANS 
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Abstract. We give a self-contained proof that the nullity of the Bezoutian matrix associ- 
ated with a pair of polynomials / and g equals the number of their common zeros counting 
multiplicities. 



IP 



With two polynomials / and g and n = maxjdeg /, deg g} we associate an n x n matrix B, 
called the Bezoutian, and a 2n x In matrix R, called the resultant. Their defining relations 
are given by © and ([5]) below, respectively. In terms of the coefficients of / and g they are 
given by ([8|) and ([!]). In this note we give a simple and self-contained proof of the equalities 



(1) 



dim ker B = dim ker R = deg gcd(/, g) , 



where gcd stands for greatest common divisor. H. Wimmer in [8] attributes this result to 
Jacobi who in 1836 showed that the singularity of what we call the Bezoutian implies the 
existence of a common factor of / and g. More contemporary proofs of (pQ) can be found in the 
recent books O Theorems 21.10 and 21.11] by H. Dym and [5j Theorem 8.30] by P. Fuhrmann. 
In the Introduction to O Chapter 21] it is shown that dim ker B > deggcd(/, g) by using the 
defining formula for B, differentiation and chains of vectors. That equality prevails is then 
proved by using these chains and the so-called Barnett identity: B = Hfg(Cj), where Hf is 
the Hankel matrix for / defined below and Cf is the companion matrix of /. In [5] the matrix 
B is expressed in terms of a matrix representation of g(Sf), where Sj is the shift operator in 
the space Xf of polynomials modulo /, relative to two suitably chosen bases in Xj. In view 
of Corollary 8.29] this formula is closely related to the Barnett identity. In this note we 
do not resort to this identity. Our approach, we think, is more direct. Of course some of the 
formulas derived below also appear in O Chapter 21] and Chapter 8]. Our proofs of these 
formulas are different. For a survey of results related to Bezoutians see [TJ Fact 4.8.6] in the 
encyclopedic book by D.S. Bernstein and for applications of Bezoutians in numerical linear 
algebra and system theory, see for example [5] and [6], respectively. 

1. Notation and basic notions. The vector space of all polynomials with coefficients in C 
and in the variable z is denoted by C[z]. Its Cartesian square is denoted by C 2 [z]. For n € N, 
C[z]< n , denotes the subspace of C[z] of all polynomials of degree strictly less than n. This space 
has dimension n. Similarly, C 2 [z] <n denotes the Cartesian square of C[z]< n . 

We use / to denote the identity matrix, Z the reverse identity and iV the nilpotent Jordan 
block: 



Z 







1 



N 



1 ••• 









For a polynomial 



/(*) = fo + flZ + ■ ■ ■ fnZ n 
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in C[z] we define two nxn matrices, one Hankel and one Toeplitz, associated with / as follows: 



Hi 



h 



fn 



in 



o 



'fo 







fn- 



fo 



Since the left-multiplication by Z reverses the rows, it turns a Hankel matrix into a Toeplitz 
and vice versa: 



ZH f 



fn 



fl 







fn 



ZTf :-- 







fo 



fo 



n-1 



TjZ, where the superscript 



As each Hankel matrix is symmetric, we have ZTf = [ZTA 
T is used to denote a matrix transpose. Consequently, 

(2) Tj = ZT f Z. 

The vector space of upper (lower) triangular Toeplitz matrices is spanned by the identity 
/ and the powers of N (N T , respectively). Therefore, the upper (lower) triangular Toeplitz 
matrices form a commutative algebra. In particular for polynomials / as above and g{z) = 
go + 9i z H 9nZ n we have 

(3) TfTg = TgTf, HfZH g = H g ZH f , 

where the last equality follows from (ZHf)(ZH g ) = (ZH g )(ZHf). 
For n £ N and z € C we denote by V n (z) the nxl column vector 

V n (z)=[l z ■■■ z n ~ l ] T . 

This notation is convenient as it provides a compact way of writing polynomials. For example, 
a polynomial a(z, w) in two variables z and w can be written as: 

n-l 

a(z,w) = ^2 a jk z j w k = V n (z) T AV n (w), 
j,k=o 

where A is the nxn coefficient matrix [%'fc]" fc _o °f cl{z,w). 

The resultant R of the polynomials / and g is the 2n x 2n matrix given as a 2 x 2 block 
matrix: 

T f ZHf 



(4) 



R 



T 9 ZH 9J 



Notice that the action of R on V-2, n (z) is particularly simple: 

(5) RV 2n {z)=[f(z)V n {z) g(z)V n (z)] T . 

Next we define the Bezoutian B of / and g. First consider the polynomial f(z)g(w) — 
f(w)g(z) in two variables. Since this polynomial vanishes for all w = z E C, there exists a 
polynomial b(z, w) in two variables such that 

f{z)g{w) — f{w)g{z) = (z — w)b(z,w) for all z, w G C. 

The Bezoutian B of f and g is the nxn coefficient matrix of b(z,w): 

(6) b(z, w) = V n (z) T BV n (w), z, w G C. 
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The null space or kernel of a matrix (or a linear transformation) A is denoted by ker A. Its 
dimension is called the nullity of A. 

2. A connection between R and B. To establish a connection between R and B we consider 
the polynomial {z n — w n )b(z, w) and we find two ways of representing its coefficient matrix. 
To find the first representation we use the standard identity 

n-1 

z n - w n = (z - w) z n ~^ j w j = {z- w)V n {z) T ZV n (w), z,weC, 

3=0 

matrix algebra and ([5]): 

(z n - w n )b{z, w) = (z- w)b(z, w)V n {z) T ZV n (w) 

= (f(z)g(w)-g(z)f(w))V n (z) T ZV n (w) 



f(z)V n (z) 
g(z)V n (z) 
f(z)V n (z) 
g(z)V n {z) 

V 2n (z) T R T 



T 



g(w)ZV n (w) 
-f(w)ZV n {w) 
z 



-Z 
z 



f(w)V n (w) 
g(w)V n {w) 



-Z 



RV 2n {w). 



The second representation involves the Bezoutian: 

(z n - w n )b(z, w) = {z n - w n )V n (z) T BV n (w) 

= (z n V n (z)) T BV n (w)-V n (z) T B(w n V n (w)) 

ro - 

V 2n (w) + V 2n (z) T 



V 2n (z) 



V 2n (z) 





B 

-B 

B 







V 2n (w) 



V 2n {w). 



These two representations of the coefficient matrix of {z n — w n )b(z, w) provide a connection 
between R and B: 









' Z~ 




ro —B~ 




(7) 




R T 




R = 












-z 0_ 




B _ 




On the other hand 


, using the definition of R, ([2]) and ([3]) 


we obtain 


R T 


" Z~ 


R = 


'ZT f Z 


ZT g Z- 




- ZT g 


H 9 ~ 






-Z 0_ 




.HfZ 


H 9 Z_ 




-ZT f 







ZT f T g -ZT g T f ZT f ZH g 



ZT g ZH f 



H f T 9 ~ H 9 T f 


H f T 9 ~ H 9 T f 



H fZH g — H g ZH 
{HfT g -H g T f ) 




•1 



BRANKO CURGUS AND AAD DIJKSMA 



Together with 0, the last equality yields 






~{HfT 9 - H g T f ) r 




"0 -B 


HfT g -H a T f 







B 



and thus 

(8) 



B 



HfT g 



H g T S = B 



3. R and B have the same nullity. Equation ([7]) indicates that there is a connection 
between ker R and ker B. An even more direct connection between ker R and ker B is obtained 
from © and (|3|) (listed in the order in which they are used) as follows: 



I 

Tf 





ZH f 



R 



0) 



T] + ZH f T g 
Tf 



ZH f 

T f ZH f + ZH f ZH g 
ZH f 

ZH g )T } {Tf + ZH g )ZHf 
I 




ZB 



Z 



I 

Tf + ZH g 
I 

T f + ZH g 



f 



ZH 





B 


o~ 




" I 










I 






ZH f _ 



If we assume that n = deg/, then Hf is invertible, yielding that the first (as well as the last) 
block matrix in (|9|) is invertible. Since the block matrix in ([9]) whose antidiagonal entries are 
/ and Z is also invertible, ([9]) implies that R and B have the same nullities: 



(10) 



dim ker R = dim ker B. 



4. The nullity of B in terms of / and g. Consider the multiplication operator 

M : C 2 [z] <n -> C[z] <2n 

defined by 



M 



fu + gv, u, v € C[z] 



<n- 



For a characterization of the null space ker M of M in terms of / and g we need the greatest 
common divisor h of / and g, its degree = deg h and factorizations / = fh, g = gh. Then 



fill 



kerM 



f 









G 



"59, v = fq, q£ C[z] <k 



The inclusion D in (jllj) is clear. To prove C, let u, v S C[z]< n and \u v] € kerM. Then 
fu + gv = 0, implies fu = —gv. Since / and g have no common zeros, the last identity yields 
that there exist polynomials p and q such that u = gp and v = fq. Substituting back to 



fu = -gv, we get fgp = -gfq. Hence p 
implies degg < k. This proves (jlip . 



-q. Since degf < n and deg/ 



n 



k, v = fq 



The standard basis for C \ 



T 




- z n-l- 




"0" 







0_ 


) * * * ) 







H- 1 


, . . . , 





while the standard basis for C[z]<2n is 



I , z, . . . . z , z 



J2n-1 
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The matrix representation for M with respect to these standard bases is R T , see ©. Therefore 
the nullity of R T is dimkerM. Since (jlip yields dimkerM = k and the nullity of R T equals 
the nullity of R, we have proved that the nullity of R is k and, by (jlOp . 

dim ker B = dim ker R = k = deg h. 

5. Final remarks. It was remarked in [HJ p. 318] that the Bezoutian of a pair of polynomials is 
defined whenever n > maxjdeg /, deg g}. We add to this that the same is true for the resultant 
and that if n > m := maxjdeg /, deg g}, then formula (H|) has to be replaced by the formula 

(12) dim ker B n = dim ker i?2n = n — m + deggcd(/, g), 

where, for example, the index n in B n indicates that B n has size n x n. Indeed, (|12j) follows 
from (pQ) and from the equalities 

dim ker B n = n — m + dim ker B m and dim ker i?2n = n — m + dim ker i?2m • 

' B o" 

The first of these two equalities holds because of ©, which implies B n = ^ m , and the 

second follows from the reasoning in Section 4 with k in (jlip replaced by n — m + k. 
Finally we note that (fl~2|) can be expressed as 

dim ker B n = dim ker R n = deg gcd(/, g) , 

where 

f(y, z) = foy n + hy n - l z + --- + f n z n and g(y, z) = g y n + g 1 y n ~ 1 z + ■■■+ g n z n 

are homogenizations (in the sense of [3 page 6-7]) of / and g, respectively. If n > m = 
maxjdeg /, deg 5}, then y = is a common zero of / and g of multiplicity n — m. Since the 
zero y = of the homogenization is commonly viewed as a "zero at infinity" of the original 
polynomial (see for example [21 4.4.3]) we can rephrase our abstract to cover the "generalized" 
Bezoutian B n : The nullity of the Bezoutian matrix B n associated with a pair of polynomials / 
and g equals the number of their common zeros including the "zero at infinity" and counting 
multiplicities. 
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